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Abstract. In the history of mechanics, there have been two points of view 
for studying mechanical systems: Newtonian and Cartesian. According the 

Descartes point of view, the motion of mechanical systems is described 
by the first-order differential equations in the N dimensional configuration 
, space Q. In this paper we develop the Cartesian approach for mechanical 

systems with constraints which are linear with respect to velocity. 
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: 1. Introduction. 

■ In " Philosophiae Naturalis Principia Mathematica" (1687), Newton consid- 

ers that movements of celestial bodies can be described by differential equations 
of the second order. To determine their trajectory, it is necessary to give the 
initial position and velocity. To reduce the equations of motion to the inves- 
tigation of a dynamics system it is necessary to double the dimension of the 
position space and to introduce the auxiliary phase space. Descartes in 1644 

proposed that the behavior of the celestial bodies be studied from another 
point of view. These ideas were stated in "Principia Philosophiae" (1644) and 
in "Discours de la metode" (1637). According to Descarte the understanding 
of cosmology starts from acceptance of the initial chaos, whose moving ele- 
ments are ordered according to certain fixed laws and form the Cosmo. He 
consider that the Universe is filled with a tenuous fluid matter (ether), which 
is constantly in a vortex motion. This motion moves the largest particle of 
matter of the vortex axis, and they subsequently form planets. Then, accord- 
ing to what Descartes wrote in his "Treatise on Light", "the material of the 
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Heaven must be rotate the planets not only about the Sun but also about their 
own centers... and this will hence form several small Heavens rotating in the 
same direction as the great Heaven." |13j. Thus the equation of motion in the 
Descartes theory must be of the first order equation in the configuration space 



Descartes gave no principles for constructing the field v for different mechanical 
systems. Hence, to determine the trajectory from Descartes's point of view it 
is necessary to give only the initial position. In the modern scientific literature 

the study of the Descarte ideas we can find in the monographic of V. V. Kozlov 
[13] in which the author affirms that "solving dynamics problem is possible 
inside the configuration space" . 

In [THl |2H |23] we developed the Cartesian approach for mechanical systems 
with three degrees of freedom and with constraints linear with respects to 
velocity. The aim of this works is to generalized the Cartesian approach for 
non-holonomic mechanical system with iV degrees of freedom and constraints 
which are linear with respect to the velocity. 

We shall present briefly the contents of the paper. 

In section 2 we prove our main results (see Theorem ll.l[ Corollary II. 3\ 
Corollary II. 4[ Corollary 11.51 below) . 

In section 3 Corollary 11.31 applied to determine Cartesian and lagrangian ap- 
proach for non-holonomic systems with three degree of freedom. We illustrate 
the obtained results to study Chapliguin-Caratheodory's sleigh and to study 
Suslov's problem for the rigid body around a fixed point. 

In section 4 we determine Cartesian and lagrangian approach in three di- 
mensional Euclidean space. 

In section 5 by applying the results of the previous section we study the 
integrability of the geodesic flow on the surface. 

In section 6 Theorem 11.11 applied to the study Gantmacher's system and 
Rattleback. 

In section 7 Corollary 11.41 applied to solve the inverse problem in dynamics. 

For simplicity we shall assume the underlying functions to be of class C°°, 
although most results remain valid under weaker hypotheses. 

It is well known that the behavior of constrained Lagrangian system 



can be described by the differential equations deduced from the D'Alembert- 
Lagrange Principle |14[ [T5] 



Q 



x = v(x, t), x G Q 




j,k=l 



k=l 



(1) 




k 
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where //i, /j, 2 , . . . , [Im are Lagrangian multipliers. 
Our main results are the following 

Theorem 1.1. Let Let Q be a smooth manifold of the dimension N with 
local coordinates x = (x 1 , x N ) and equipped by the Riemann metric G = 
(Gkj(x)) = (Gkj) and let 



N 



a = ( v(x) , dx) = ^kj 
j,k=i 

be the 1-form associated to the vector field 



G kj x j v k , 



(2) 



v 



T 





^i{d 2 ) ■ 


■ ^i(d N ) 







^m(<9 2 ) • 









^M+l(9 2 ) • 


■ ^m+i(9n) 


Am+i 






■ M N (d N ) 


Xn 




d 2 


d N 






d 

where T = Q 1 A Q 2 ... A Q N (di, d 2 , djv), <9 fc = 7^, Aj = Aj(x) /or j = 

M + 1, . . . N are arbitrary functions and Q k for k = 1, ... N are 1-forms on Q 
which we assume that satisfies the following conditions 

N 

(i) Qj for j — 1, . . . , M are a given 1-forms: Qj = ctjkdx k where = 

k=i 

OLjk{x) are functions on Q, 

(ii) Qk for k — M + l, . . . , N are arbitrary 1-forms which we choose in such 
a way that T / 0. 

(iii) The 2-form da admits the development 

1 N 

da = - ^2 ajk{x)Qj A Q k 
j,k=i 

where A = (a^) is a skew symmetric matrix such that 

H = M T AM 

where H = (da(dj, dk)) and M = (Qj(dk)) are N x N matrix. 

(iv) The contraction of the 2-form da along the vector field v is such that 

N 

i v da = Aj(x)Qj, 

3=1 

where A = (A 1; . . . , Ajv) t which we can be calculated as follows 
(3) A = AX = M-V, 

where X = (Ai, . . . , Xn) and r = (i v da{d\) , . . . , L v da{dN)) 
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Then the first order differential system on Q 

(4) x= v(x), under the conditions, A-^x) = 0, 

for j = M + 1, . . . N, is invariant relationship of the second order differential 
system 

Comparing equations ([I]) and (jSJ) we deduce that the latter can be inter- 
preted as the equations describing the behavior of non-holonomic mechani- 
cal systems under the action of active forces with potential U = -||v|| 2 + 

Uq, Uq — const and under the action of the reactive forces with the compo- 
nents 

(MM M \ 

j=l j=l i=l ) 

N 

generated by the constraints fij(x) = ctjkX k = 0, for j = 1, 2, ...M. 

fc=i 

Of interest is that the equations Aj = for j = M + 1, . . . , N or, which is 
the same 

N 

(6) Aj = OjfcAfc = 0, ajk = —akj 

k=i 

for j = M + 1, . . . , N, represent a system of partial differential equations of 
first order with respect to the functions for k = 1, . . . , N. 

Definition 

We call the vector field v which generated system @ Cartesian vector field. 
The vector field v we say Cartesian equivalent if there exist a nonzero function 
k on Q such that kv is Cartesian vector field. 

Definition 

Studying the behavior of nonholonomic mechanical systems with constraints 
linear with respect to the velocity using the equations (JT]), @ and (jSJ) we called 
Classical, Decartes and Lagrangian approach respectively. 

Conjecture 1.2. There are solutions of the equations ([6]) which generate a 
Cartesian (or Cartesian equivalent) vector field which completely describe the 
behavior of the study constrained Lagrangian system. 

This conjecture supports the following facts. 

First, in view of Theorer dl.il the solutions of @ are solutions of (jSJ), which 
is closely linked to the system (JT]). Second, the solutions of the equations (JTJ 
depend on the 2A^ — M initial conditions. The solutions of (jl]) depend on 
N initial conditions and N — M functions which are solutions of the linear 
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partial differential equations ((6]), therefore the solutions of the equations (EJ) 
also depend on N — M arbitrary constants. Finally, also contribute to the 
strengthening of the conjecture the large number of applications given below. 

Corollary 1.3. Let us suppose that in Theorem li.il manifold Q is three di- 
mensional smooth manifold with local coordinates x = (x, y, z) and the given 
1-form is Q = a\dx + a 2 dy + a 3 dz = (a, dx) where aj = aj(x, y, z) for j = 1, 2, 3 
are functions on Q. 

Denoting by j x ] the vector product in R 3 and by rotv, a, w the fol- 
lowing vectors fields 

rotv = J- — ((d y p 3 - d z p 2 ), (d z pi - d x p 3 ), (d x p 2 - d y pi)) T , 



VdetG 
1 



w 



T 



((A 2 ^ 3 - A 3 ft 2 )(3r), (A 2 ^ 3 - x 3 n 2 )(d y ), (\ 2 n 3 - x 3 n 2 )(d z )y 



where pk = GkjV-* , for k — 1, 2, 3, then differential system (j4|) and (jSJ) take 



the form respectively 



(7) 



x = [a x w] = v(x), (a, rot[a x to]) = 0. 



dd,T 8T d 2^ v \ 



w, rot[a x to]) a 



dt dx dx dx 

— £l 2 A Q 3 (v, rotv) a 



d-\\v\ 



dx 



d 



d 



where — = (d x i , . . . , d ±N ) , — = (d x i, . . . , d, 
dx dx 

Corollary 1.4. Let us introduce the notation 

d/i(0O ... dh(d N ) 



... df N -x(d N ) 
di ... d N 



{fi, fa, /jV-l) *} 



and suppose that in Theorem \1.1\ the given independent 1-form are such that 
Qj = dfj(x) for j = 1,2, N — 1, and the arbitrary 1-form it is also exact, 
i.e. tt N = df N , £l N (v) = X N , and such that T = {/ ls fa, f N _ u f N } ^ 0. 
Then the equations (jl]) and <§5§ take the form respectively 

{fi, fa, • • • , /at-i, *} 



x 



-A 



N 



(9) 



{fi, fa, ■ ■ ■ , /at-i, /at} 

1 



A{/i, fa, ■ /at-i, *}, 



d dT dT _ d ^ 



dt dx k dx k 



dx k 



2 N-l 

- + A ^ a Njdfj(d k ), 

3=1 
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where Xn = Xn(x) is an arbitrary function and a^j = aNj(x) for j = 1, ... N— 
1 are element of the skew symmetric matrix A. 

Corollary 1.5. Differential equations (jSJ) are Lagrangian with Lagrangian 
function L = —\\x— v(x)\\ 2 , where v is Cartesian vector field. 



The proofs of Theorem 11.11 , Corollary 11.31 Corollary 11.41 and Corollary 11.51 
are given in section [2J 



2. Proof of the main results 



Proof of Theorem \l.ll Firstly we shall introduce the following notation and 
concepts. Let £ G (Q) be the Lie algebra of vector fields on Q and let V be 
the connection 

V:£(Q) x£(Q) •—►£((}), 

(U, V) I ► V u v , 

which is R lineal with respect to v and C°° lineal with respect to u and is 
compatible with metric G, i.e. V u G(v, w) = 0, Vu, v,w 6 £(Q). 
Let v e £(Q) be a vector field: 



T 



n 1 (d 1 ) 




... n^dx) 


Ai 








A 2 






■ ■ ■ Qm+i(9n) 


Am+i 












d 2 


d N 






where T = f2i A Q 2 - - A ^^(81,82, On) 7^ 0. 

The functions Xj for j = 1, . . . N, are arbitrary functions on Q such that 
fy.(v) = Xj, for j = 1,...N. 

Let a be the 1-form associated with the vector field v, i.e. 

N N 

a = (v(x),dx) = GjkV{x)dx k = '^^p k dx k , 

j,k=l k=l 

then, in view of the condition T ^ the 2-form da admits the development 
1 N 

da = — a,jk(x)£lj Afifc, where A = (ttjk) is a N x N skew-symmetric matrix 

j,k=i 
such that 

a jk = (-ly+k-i—da A Q 1 A ... A Q k .... A Q r ... A fi^di, 9 2 , 9jv), 



fL, fl k , means that these elements are omitted. 
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In view of the relations 

. N N N 

Lvdcr = - ^2 (%fc^j(v) - a fe jOj(v)) Q k = ^2 a jk Qj(v)tt k ,= a jk \jQ k 

k,j=l k,j=l k,j=l 

we obtain that the contraction of da along v is 

N 

(10) tvd5- = 5^A.A', 

3=1 

where A = ^Ai, A2, ...,Anj . It is easy to check that this vector can be cal- 
culated as follows 



;il) A = AX = M r, 



where M = (Qj(d k )) and A = (Ai, . . . , Xn) T , f = (i^da(di) , . . . , L-vda(dN)) T ■ 
Now we prove that the differential system 

(12) x = v(x), xeQ, 

is invariant relationship of the Lagrangian equations with Lagrangian function 

1 1 N 

L= -||x- v|| 2 = - G kj {V- V){x k -v k ). 

j,k=i 

Indeed after covariant differentiation we obtain V x (x — v) = 0, or, what is 
the same 

V X G (x - v) = V x = 0, 

hence , by considering that 

ddT dT 1 
Vx(Gx)= dtd±~^ T= 2 I|X|1 

N N 

3=1 3=1 
V N N 

^x 3 (djp k - d k pj) + J2v J V d Jj + ^2{x J - v 3 ) V dk pj 
3=1 3=1 3=1 

N 



" iiv "'^P 



dt dx k dx k 

3=1 



N 



where p k = S ^G k jV 3 \ V = (x,v) — -||v|| 2 , and ( , ) is the scalar product. 
3=1 

d dL dL 

Then along the solutions of ( 1T2"1) we give -tttt ~ ttt = 0, for A; = 1, 2, ... AT. 

dtox h ax k 
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It is easy to show that these Lagrangian equations admits the representation 

d dT dT 



where u 



dt dx k dx k 

|v|| 2 + Lvda. In view of ( TTTJT) and ( |T2l we finally deduce the 



differential equations 

d dT 

(13) 



dT 



d-\\w\\ 2 

2 M II 



A' 



3=1 



dt dx k dx k dx k 
If we determine the vector field v and functions Ai, A 2 , . . . as follows 
(14) v = v| Ai =a 2 =...=a m =o, Aj = A i | Al= A 2 =...=A M =o, 3 



1,2,. ..N 



and require that Aj = 0, j = M + 1,...,N, then we obtain the Cartesian 
vector field and differential system (TIB"]) coincide with system (jSJ). In short 
Theorem 11.11 is proved. □ 

Proof of Corollary \1.3[ The vector field ([2]) in this case takes the form 



Tv 



ni(d x ) ni(dy) n^) o 
n 2 (d x ) n 2 (d y ) n 2 (d z ) a 2 
n 3 (d x ) n 3 (d y ) n 3 (d z ) a 3 



9a 



On 



o 





n x {d x ) 






n 3 (d x ) 






d x 








nx(9 v ) 


A 3 








d x 








Oi 


a 2 a 3 






u>i 


W 2 ^3 

<9 y d z 



d z 



thus v(x) = [a x w]. 

On the other hand considering that 

da = (rotv)a;dz A dy + (rotv) y <ix A + (rotv)^dy A dx 

where (rotv) x = cfx(rotv), (rotv) y = dy (rotv), (rotv) 2 = dz(rotv), we get 

t v da = ([v x rotv], dx) = AiOi + A 2 fi 2 + A 3 fi 3 
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which lead to 

TAi = Q 2 A n 3 (v, rotv) = A 2 fi 3 (rotv) - A 3 fi 2 (rotv) = 

= (w, rotv) , 

TA 2 =fi 3 A fii(v, rotv) = A 3 fii(rotv) - Aifi 3 (rotv), 

= A 3 ^i(rotv) = A 3 (a, rotv) , 

TA 3 = -fi 2 A Qi(v, rotv) = -A 2 Oi(rotv) + Ai^ 2 (rotv), 

= — A 2 f2i(rotv) = — A 2 (a, rotv) , 

here we put fii(v) = Ai = 0. 

From the conditions A 2 = A 3 = we obtain 

(a, rotv) = (a, rot [a x w]) = 0, 

hence we easily deduce differential system ([71) and (jHJ). 
In short Corollary 11.31 is proved. 



□ 



Proof of Corollary TJj. In this case we obtain that the vector field v takes the 
form 

dh{d x ) ... dh{d N ) 



1 

T 



df N -i{di) 
df N (di) 



df N -i{d N ) 
dfN{d N ) X N 
d N 



Ajy 



dh{dx) 



dfi(d_ 



N 



df N _i(d N ) 
d N 



-A 



JV" 



{h, h ■ ■ ■ , /at-i, *} 



A{A, h . . -Jn-i, *}■ 



{fl, f'2 ■ ■ ■ , /jv-1) /at} 
On the other hand by considering that Xj = for j = 1,2, ... N — 1, from 
03]) follows that Aj = a^jXiy, j = 1, . . . , N — 1, A N = o^nXn = 0. 

Clearly the last equation is satisfied identically in view that (a^) are ele- 
ments of the skew-symmetric matrix A. 

Therefore we easily deduced the differential equations generated from Carte- 
sian and Lagrangian approach. 

Thus Corollary 11.41 has been proven. □ 

Proof of Corollary 11.51 Follows from the proof of Theorem 11.11 by considering 

can □ 



3. DECARTES APPROACH FOR NON-HOLONOMIC SYSTEM WITH THREE 
DEGREE OF FREEDOM AND ONE CONSTRAINTS . 



In this section we apply the corollary II. 3l to study the Chapliguin-Catatheodory 
sleight and Suslov's problem for the rigid body around a fixed point. 
Chapliguin-Carathodory's sleigh 

We shall now analyze one of the classical nonholonomic systems Chapliguin- 
Carathodory's sleigh (which we call a sleigh). [TH] 
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The idealized sleigh is a body that has three points of contact with the plane. 
Two of them slide freely but the third, A, behaves like a knife edge subjected 
to a constraining force R which does not allow transversal velocity. More 
precisely, let yoz be an inertial frame and £ Ar] a frame moving with the sleigh. 
Take as generalized coordinates the Decartes coordinates of the center of mass 
C of the sleigh and the angle x between the y and the £ axis. The reaction 
force R against the runners is exerted laterally at the point of application A 
in such a way that the r] component of the velocity is zero. Hence, one has the 
constrained system with the configuration space Q = S 1 x IR 2 , with the kinetic 
energy T = — (y 2 + z 2 ) + y£ 2 , and with the constraint ex + smxy — cosxi = 0, 
where m is the mass of the system and I c is the moment of inertia about a 
vertical axis through C and e = \ AC\. 

Observe that the "javelin" (or arrow or Chapliguin's skate) is a particular 
case of a sleigh and can be obtained when e = 0. 

To apply the Decartes approach for this system, first we introduce the 1-form 
Qj, for j = 1,2, 3 in such a way that the determinant T ^ 0. This condition 
holds in particular if 

Qi = edx + sin xdy — cos xdz, Q 2 = sin xdz + cos xdy, f2 3 = dx, 

so T = Hi A 2 A ^3(<9ie, d y , d z ) = 1. 

The Descartes approach produce the differential equations (see formula ([7])) 



here Xj = A 3 (a;, y, z, e) for j = 2, 3 are solutions of the equation 

(16) sinx(J<9 2 A3 + em<9 y A 2 ) + cosx(Jd y \ 3 — em<9 z A 2 ) — m(<9 x A 2 — eA 3 ) = 0, 

where J = J c + e 2 m, \ |v| | 2 = JA 2 + mA 2 . 

Now we show that there are solutions of (fT5l) and ([16]) fully describes the 
inertial movements of the sleigh. 

Corollary 3.1. All the inertial trajectories of Chapliguin -Carathodory's sleigh 
can be obtained from Cartesian approach. 

Proof. Let us suppose that Xj = Xj(x, e) for j = 2, 3. Clearly that in this case 
(TT6"j) takes the form d x X 2 — eA 3 = and all paths of the equation (fTS]) can be 
obtained from the formula 



[24j 



(15) 



x = X 



y = A 2 cos x — eX 3 sin x, z = A 2 sin x + eA 3 cos x, 
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On the other hand, the inertial motions of the sleigh can be obtained from 
the equations deduce from the classical approach 

Jcx = e/i, my = fismx, mz — — /icosx, ex + sinxy — cosxz = 0. 

Hence, after straightforward calculations we get the first order ordinary differ- 
ential equations 

x = qC cos 9, y = C (sin 9 cos x—qe cos 9 sin x) , z = C (sin 9 sin x+qe cos 9 cos x) , 

where 9 = qex + C, q = 

This system can be obtained from (1151) if we choose 

A 2 = Csin#, A 3 = Cq cos 9. 

Is clear that in this case 

2||v|| 2 = JXl(x,e) + m\l(x,e) = mC 2 , 
therefore the sleigh moves by inertia. So the corollary is proved. □ 

Now we study Chapliguin's skate. Cartesian approach in this case produce 
the differential equations, which can be obtained from (fl5l) and (fl6|) by putting 
e = 0. 

x= \ 3 (x,y,z,0), y = \ 2 (x,y, z,0) cosx, z = X 2 (x, y, z, 0) sinx, 

J(sinxd z \ 3 + cosxdyXs) — md x \2 = 

Corollary 3.2. All the trajectories of Chapliguin's skate with the initial con- 
dition x(to) = Co 7^ and under the action of the potential field of force with 
potential function U = mgy can be obtained from the Cartesian approach. 

Proof. In fact, for the case when e = the classical approach for Chapliguin- 
Carathodory's sleigh gives (Chapliguin's skate) the following equations of mo- 
tion 

x = 0, y = g + /isinx, z = — /icosx, sin xy — cos xz = 

Hence, integrating we deduce the following differential system of first order 
(see for instance ISM) 



x = C 0, y = (— hCijcosx, z — (— — hCijsmx. 

Co Cq 

Let b be the vector field associated with this differential system, i.e. 

gsmx gsmx . 
Co, {— hCijcosx, (— — hCijsmx 

^0 Co 

9 cosx 

Whereas (a, rotb) = — — — ^ o. 

Co 
C 

Denoting by k = — — — , we easily obtain that (a, rot(« b)) = 0, 

gsmxGo + Gi) 

so b is Cartesian equivalent vector field. 



12 
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The vector field Kb can be obtained from (IT71) by choosing A3 = k and 

A2 = Co- 
Summarizing the corollary is proved. □ 

The rigid body around a fixed point in the Suslov case 

in this section we study one classical problem of non-holonomic dynamics for- 
mulated by Suslov [251 [H]. We consider the rotational motion of a rigid body 
around a fixed point and subject to the non-holonomic constraints (a, u) = 
where 00 = (co>i, Co>2, ^3) is a body angular velocity and a is a constant vector. 
Suppose the body rotates in an force field with potential Z7(7i, 72, 73)- Apply- 
ing the method of Lagrange multipliers we write the equations of motion in 
the form 

Iuj = [Iu x u] + [7 x f£] + /ia, 7 = [7 x u] 

(18) 

(a,w) = 

Where 7 = 71,72,73) = (sin z sin x, sin z cos x, 73 = cos z)„ I is the inertial 
tensor of the body, \x is the Lagrange multiplier which can be expressed as a 
function of u and 7 as follows 

(a, [Iu, qj + [7, g]) 
(a, /-ia) 

It is well-known the following result [2] 
Proposition 3.3. If a is an eigenvector of operator I, i.e. 
(19) la = na, 

then the phase flow of system (|18[) preserves the " standard" measure in M 6 = 
M 3 {^} x R 3 { 7 }. 

G.K. Suslov has considered a particular case when the body is not under 
action of exterior forces: {7 = 0. If f|T9|) holds then the equations (|18p have 
the additional first integral K4 = (Iu,Iu). E.I. Khar lamova in [11] study the 
case when the body rotates in the homogenous force field with the poten- 
tial U = (b, 7) where b is an orthogonal to a vector. Under these condi- 
tions the equation of motion have the first integral K± = (Iu,b). V.V. Ko- 
zlov in [13] consider an opposite case, when b = ea, e 7^ 0. The integrabil- 
ity problem in this case was study in particular in [TH [17] . The case when 
U = edet/(/ _1 7, 7) the system (|T8l) have the Clebsch-Tisseran first integral 

K 4 = - -edet/OT-S, 7), M- 

From now on, we suppose that equality ( TT9l) is fulfilled. We assume that 
vector a coincides with one of the principal axes and without loss of generality 
we can choose it as the third axis, i.e., a = (0, 0, 1) (see for more details |14j ) 
Equations of motion have the following form 
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iiWi = 73^ U - 7 2 <9 73 U, I 2 u 2 = 71^73 u - 73^ 71 U, 
(20) { (/ 1 -/ 2 )w 1 w 2 + 7 2 9 7l ^-7ia 72 f/ + ^ = 0, 

7i = -73^2, 72 = 73^1, 73 = 7i^2 - 72^1, 

where Ii, I 2 are the principal moments of inertia of the body with respect to 
the first and the second axis. We observe that the above mentioned choice of a 
guarantees that the phase flow of system (120]) preserves the standard measure 
in M 5 {u;i, u 2 , 7}. 

In [20] we prove the following result 



Theorem 3.4. Let us suppose that the body rotates within the force field defined 
by the potential 

(21) U=-± T {h& + h£)-h, 

where h is a constant and fii, \x 2 are the solutions of the following first order 
partial differential equation 

22 73- — - 72 _ + 7l _ = 0, 

^72 9ti ^73 £73 

then the following statements hold: 

(I) The equations (1201) have the first integrals 

(23) hui — H2 = 0, I 2 u 2 + hi = 0, 

consequently, they are integrable by quadratures. In particular 

55(7!, 72, K 2 ) 2 2 

huJi= ~ + *2(7i +73^2,72) 

C72 



(24) 



+7 1 fi(7 1 2 + 7 2 i ,K 2 ,7 3 ), 

95(7i, 72,-^2) , Tf , 2 . 2 v v 

^2= 5 *lW + 73^2,7l) 

071 



-72^(7^ + 7 2 ,if 2 , 73), 

are constants on the solutions of (1201 . where K 2 = 7i + 7 2 + 73 awe? 
5,^1, ^2, ^ are arbitrary smooth functions . 
(II) The Suslov's, Kharlamova-Zabelina's, Kozlov's , Clebsch-Tisseran 's and 

Tisseran-Okunova's first integrals can be obtained from (1241) . 
(Ill) The dependence 7 = j(t) we determine by quadratures of the Poisson 
equations which in this case take the form: 

fnK \ • 73^i . 73/^2 . 7i/*i 

(25) 7l = ——, 72 = —j—, 73 = = r—. 

1 2 H 1 2 H 

It is interesting to note that the proof of the Theorem 13.41 first was obtained 
using Cartesian approach (for more details see [21]) which we proposed below 
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Let us suppose that manifold Q is the special orthogonal group of rotations 
of E 3 , i.e. Q = 50(3), with the Riemann metric 

(I 3 I 3 cos z 

I 3 cos z (Ji sin 2 x + I 2 cos 2 x) sin 2 z + 1 3 cos 2 z (h — I 2 ) sin x cos x sin z 
(ii — J2) sin x cos x sin 2; ii cos 2 x + I 2 sin 2 x 

det G = hhh sin 2 2. 

In this case we have that the constraint is u 3 = x + cos zy = 0. 

Hence a = (1, cos z, 0). By choosing the 1-form Qj for j — 1, 2, 3 as follow 

f2i = rfx + cos z<i?/, ^2 = ^2/, ^3 = dz. 

Consequently T = dVL\ A dQ 2 A dfi 3 (9 x , 9y, d z ) = 1, and 

v = A 2 (cos zd x - d y ) - X 3 d z 

Thus we obtain that 
Pi = 0, 

P2 = (^3 _ h + (A _ ^2) cos 2 x) cos 2; sin 2 z\ 2 + (ii — ^2) cos x sin X sin z\ 3 , 
P3 = (^2 sin 2 x + Ii cos 2 x) A3 + (I 2 — h) sin x cos x sin 2 cos z\ 2 

The differential equations ([7]) in this cases take the form respectively 



(26) x = cosz\ 2 , y = -\ 2 , z = -A 3 
and 

(27) (a, rotv) = <9 2 p 2 - d y p 3 + cos = 

After the change 71 = sin z sin x, 72 = sin z cos x , 7 3 = cosz the system 
( 12"B"|) and condition (12T|) can be written as follow 

1 1 1 

(28) 71 = — ^173, 72 = ^-^273, 73 = — p— - (ii/ii7i + 

J-2 *1 -'1J2 



(29) smz 73 -— ) - 7 2T ^ h 7it^— ) - cosxd y n 2 - smxdy^ = 

£72 5 7i ^7s ^73 

where /z 2 = — Ii(cosxA 3 + sinxA 2 ), 1^1 = I 2 {— sinxA 3 + cosxA 2 ). 

Clearly if ), for j = 1,2, then the equation (1291) coincide with 

equation ([22]) and (EH]) coincide with ([25]). 
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4. Cartesian vector field in three dimensional Euclidean space 

Let E 3 be the three dimensional Euclidian space with Cartesian coordinates 
x = (x 1 ,x 2 ,x 3 ). 

We consider a particle with Lagrangian function 

L = i||x|| 2 -t/(x), 

and constraint (a, x) = 0, where (, ) denotes the scalar product in E 3 , x = 
(xi, X2, is) and a = a(x) = (ai(x) ,0,2(2:), a 3 (x)) is a smooth vector field in 
E 3 . Below we shall use the following notation <9 X / = (d xi f, d X2 f, d X3 f) T . 

The equations of motion in particular for constrained particle in R 3 , can be 
deduced from the d Alembert-Lagrange Principle and are such that 

x = dyU + /ia, (a, x) = 0, 

where \i is the Lagrangian multiplier. 

Cartesian and Lagrangian approach produces the following differential equa- 
tions respectively (see formula (J7|) and OH])) 

x= [axw], (a, rot [a x w]) = 0, 
(3°) .. n ( 1 , 



<9 X Q||[a x w] 1 1 2 ^ + (rot [a X w], w)a, 



Example 

c 

Suppose that a = f x , w = — , 

c 2 

where / = r + (x, b) = c 2 and b and c are constants vector field such that 

(x,c) = 0, (b,c) = 0, ||c|| 2 = c 2 , r = ||x||. 

Then Cartesian and Lagrangian approach generate the following differential 
equations respectively 

± = [jx X C] = X 

c 

Indeed in view of the relation rot[/ x x c] = — , we get 

r 

(a,rot[a x w]) = (/ x , = 
thus the given vector field is Cartesian. 

From the relations (rot [a x w],w) = — j(c, rot[/ x x c]) = -, we 

obtain that the second order differential system ( 130]) in this case takes the 
form 

1 ( o J, d*f\ f x x 



x 



4 dJ^) - 
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Below we study the particular case when the vector field a satisfies the 
equation 

(31) rot(a) = u(x)a., 

where v is certain function. 
Optical-mechanical analogy 

From the standpoint of geometric optics, propagation of light in E 3 can be 
represented as a flow of particles. Trajectories of particle are called rays. 
It is known [TJ] that the vector field K of an arbitrary system of rays an a 
homogeneous optical medium satisfies the relation K x rotK = 0. System of 
rays such that rotK ^ are called Kummer systems 

Proposition 4.1. Let a be the Kummer vector field i.e., satisfies the partial 
differential equation ( 13 ip with v ^ 0. Then Cartesian and Lagrangian approach 
generate the following differential equations respectively 



(32) x = V- 7a = jT—r^la x rotW\ = v(x) 

\\ a \\ 

and 

(33) x=d x (h [ \[ax HII 2 ) +{&,[ ax «>])) <*» 

where 7 = - ' ^ , V and R are the vector fields: 
\\ a \\ 

(34) rotV=uV+[ax V\, R=V-d x y, 
where V is an arbitrary smooth vector field. 

Proof. Indeed, taking the well-known relations 

div[A x B] = (A,rotB) - (rotA,B). 

into account, we obtain 

div[a x [a x w]] = (a, rot [a x w]) — (rota, [a x w]) 

Hence in view of this identity, and considering that (a, rot [a x w]) = and 
( 13~T|) we obtain that [a x [a x w]] is a solenoidal vector field, so 

(35) [ax [ax w]] = (a, a)w — (a, w)a = rotW, 

where W is a smooth vector field such that (a, rotW) = 0, consequently 
rotW = [a x V], for a smooth vector field V which must satisfy the partial 
differential equation 

div(rotW) = div([a x V]) 

= (a, rotV) - (V, rota) = (a, rotV - uV) = 0, 



hence we obtain the representation ([34 
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In view of relation ( 1351) we get the following representation for w 

(a, w) [a x V] (a, w) rotW 
w == — — rrr-a H — t: — r— — = — — rrr^-a H — - — — 
\\a.\\* ll a ll ll a ll ll a ll 

From the relation v = [a x w] we get formula f l32l . 
On the other hand, after some calculations we obtain 
rotv = rotV — [<9 X 7 x a] — 7rota 
= v (V — 7a) + [V — 9 X 7 x a] 

or, what is the same 

(36) rotv = uv+ [R x a]. 

Inserting (1361) in the relation (w, rotv) we obtain 

(w, rotv) = (w, w + [R x a]) 

= (w, [R x a]) = (— [ w x a], R) 
= ( [ a x w] , R) . 

Hence from (130|) we obtain the second order differential equations (|33|) . □ 
Corollary 4.2. Let us suppose that Cartesian vector field is such that 

(37) rotv = i/v, v 7^ 0, 

then Lagrangian approach generated the differential system 

1,, 



x= d x 

which describe the motion of a material point of unit mass in the potential field 
with the force function — ||v|| 2 . 

Proof. From ( 1371) and ( 1361) follows that [R x a] = 0, based on this relation and 
considering that 

(R, [a x w])) = (w, [R x a])) 

we deduce that (R, [a x w])) = 0, thus, in view of (1331) . we obtain the proof 
of the proposition. 

This results coincide with theorem J.Bernoulli (1696.) " Light rays in an 
isotropic optical medium with the refraction index n(x) coincide with the 
trajectories of a material point in a potential field with the force function 

Example 

Consider a particle with Lagrangian function L = ^i^ 2 + V 2 + i; 2 ) and 
constraint o>2 = sin x y — cos x z — 0, 

The vector field a = (0, sinx, — cosx), satisfies the equation rota = a. 
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Classical approach generates the differential equations 
x = 0, y — fisinx, z= — ficosx. 
Integrating these equations we obtain 

(38) x — Ci, y = C 2 cosx, i = C 2 sina;, 

where C\ and C 2 are arbitrary constants. The vector field associated to this 
system is Cartesian type. 

Cartesian approach generate the differential system 

x = X, y = g cosx, z = g siiax 

where A = u>3 sin x + w 2 cosx, g = —W\, and W\, w 2 , u>3 are components of the 
arbitrary vector field w. 

From the condition (a, rot [a x w]) = follows that the functions A and g 
are solutions of the linear partial differential equation 

(39) d z X sin x + d y X cosx — d x g = 0. 
In particular the vector field 

(40) v= (X(x), g(y,z) cosx, g(y, z) sin x) , 

satisfies (13"9"j) . where A = X(x) and g = g{y,z) are arbitrary real functions, . 

Solving equations A (a?) = w 3 sin x + w 2 cos x, g(y, z) = —W\, with respect to 
the components of the vector w we obtain 

w = (—g(y, z), X(x) cosx — ifi(x, y, z) sinx, X(x) sinx + ip(x, y, z) cosx) 

where -0 = ip(x, y, z) is an arbitrary function on Q. 

It is easy to check that 

rotv = v + u, 

where u = {d z gcosx — d y gsmx — X, 0, 0) T , is a vector orthogonal to vector a. 
Hence we obtain the relations 



[Rxa] = u, (w, [Rxa]) = g (d y gsinx + A — d z gcosx) , ||[wxa]|| 2 = g 2 +X 2 . 
Lagrangian approach generates the equations 

x = 9 X {^{q 2 + ^ 2 ) S j ~ Q (d z gcosx — d y gsmx — A) a 

The vector field ( I38jl can be obtained from (HU1) if we choose X — C\, g — 
C*2. For these parameter values Lagrangian approach generates the differential 
equations x = CiC 2 a. 

If we choose g and A: (d z gcosx — d y gsinx — A = Then rotv = v and as a 
consequence Lagrangian approach generate the equation x = 5 X f - (g 2 + A 2 )^ . 

Below we study the case when the constraints are integrable. 
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Corollary 4.3. Let V and a are the vector field such that 

(41) rota=0, V=d f Gd x <S> 

where G = G(f, $) and $ are an arbitrary smooth functions, then Cartesian 
and Lagrangian approach for a particle in E 3 which is constrained to move on 
the surface f = f(x) = c generate the following differential equations 

(42) x=d f G(d x $- 1 d x f) = v(x), 
and 

f d G 2 

^ x=d x [^-\\d x f xd x ®\ 

-d f G {d} f Gd x § - d xl , d x $ - !d x f) d x f, 

where g = \ \d x f\\ 2 , and 7 = - — — - — , 

9 

7/7 = then the equations ( 1421) and (1131 take the form respectively 

(44) x = d f Gd x <f>, ■x=d x (^\\dM 2 ) -(G ff G\\dM 2 )dJ 

Proof. Suppose that (HIT) hold, then 

a = d x f, rotV = [d*f x (d 2 ff Gd^) = [d*f x V], 
, hence the vector fields W, V and R admit the representation 

W = G<9 X $, V = d) f Gd^ + vdxf, R = d 2 ff Gd^ + vd x f - d x7 , 
where v is an arbitrary function, hence we easily obtain (I43p . 

If 7 = then (R, d f G (<9 X $ - 78*/)) = d f Gd 2 ff G\\d x <5>\\ 2 , thus formula (l44j) 
follow . □ 

5. INTEGRABILITY OF THE GEODESIC FLOW ON THE SURFACE. 

It is well known that the differential equations x = /i<9 x /, where \x is La- 
grangian multiplier, determine the geodesic flows on the surface / = c, and 
admits the energy integral 

||x|| 2 = 2M/). 

If there is an additional first integral, functionally independent with the energy 
integral , then the geodesic flow is integrable. 

Proposition 5.1. Let us suppose that ( |4"T|) holds. Then Lagrangian geodesic 
flow of the constrained particle on the surface f = c, g = ||<9 X /|| 2 > is 
integrable if there exist a solution of the following non-lineal partial differential 
equation respectively 

(45) G 2 f (f^)\\id x fxdM\ 2 = 2Hf)g, if (d x f,d x ®)^o, 

and 

(46) G 2 f {f,$)\\dM 2 = 2h(f), if {d x f, d x $) = 0, d x $^KX. 
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Proof. Indeed, if ( 1451) holds then the system ( 1431) takes the form 

(47) x = (d f h(f) - d f G {8P ff Gd x $ - d xl , a x $ - jd x f)) d x f, 

which determine the geodesic Lagrangian flow and admits the following com- 
plementary first integral 



4||[d x $xx]|| 2 
1 (d x f, d x $) 2 Ul ' 



which it is easy to obtain from (I42p . Clearly this first integral is independent 
of energy integral. 

If (d x f, <9 X $) = then Cartesian and Lagrangian approach generated the 
differential equations respectively 

x = G f (f, $ )<9 X <&, x — (hf — G f (f, $ )Gff(f, $ )||9 X $| | 2 )9 X /, 

The condition (145 p under the given condition of orthogonality takes the form 



C 2 , 



The complementary first integral is 

||$ x || 2 ||fx x xl 



H[xxa x $]|| 2 

which we can obtain from ( 1421) . by considering that 7 = 0. This complete the 
proof of the proposition. □ 

Now we apply the above results. 
Now we consider the surface 

(48) f{x)=c, (x,a K /) = m/, c^O 

which we call homogeneous surface of degree m. 

We are interested in studying the integrability of Lagrangian geodesic flow 
on the homogenous surface. 

Euler's formula shows that c = is the unique critical value of /, hence for 
c 7^ the function g = ||<9 X /|| 2 > 0, on the surface f(x) = c. 

Taking into account formula (|48p it follows that 

(x, d x g) = 2(m - l)g. 

Below we use the notation 

d x F d y F d z F 



{F, G, H} 



d x G d y G d z G 
d x H d y H d z H 



Clearly, if F, G, H are independent functions then {F, G, H} ^ 0. 

The integrability of the geodesic flow on the homogeneous surface we study 
in the following two cases 

(49) {/,4,r 2 } = 0, {f,g,r 2 }^0, 

where r 2 = x 2 + y 2 + z 2 . 
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We analyze the first case. We study only the particular subcase when the 
homogeneous surface satisfies the condition 

(50) g = g(f,r). 

Hence, in view of (jUJ) we obtain mfdfg + rd r g = 2(m — l)g. 

We assume that function $ is such that $ x = x, thus the differential equation 
generated by the Cartesian and Lagrangian approach are respectively 

(51) ±=£V-m/<W), S=!=*fW>«J, 

g r z g z 

Proposition 5.2. Lagrangian geodesic flow on the homogeneous surface under 
the assumption f loUj) is integrable 

Proof. Let us suppose that <9 X $ = x, then (<9 X $, d^f) = mf ^ 0. On the other 
hand if we choose G(f, r) as 

C 2 (/;r) _ Mf)9(M 

' g{fi r ) r2 ~ m 2 f 2 ' 
then we obtain that (1451) holds. Thus there exist an additional first integral 

s(/,r)||[xxx]|| 2 = m 2 A(Z) 

□ 

Example 

Lagrangian geodesic flow on the homogeneous surface of degree one 

f(x) = r + (b, x) = c, c ^ 
is integrable, where b is a constant vector field. 

2/ 

Proof. In this case we have g = h llbll 2 — 1 = g(f,r). 

r 

The complementary first integrals is 



(^ + ||b|| 2 -l)||[xxx]|| 2 = 2A(/). 



r 



□ 



We have studied the case in which {f,g, r 2 } = 0. Now we study the case 
in which the functions /, g, r 2 are independent, i.e., {/, g,r 2 } ^ 0. Hence, we 
obtain that 

(52) x = x(r,f,g), y = y(r,f,g), z = z(r,f,g). 

To establish the integrability or non-integrability of the Lagrangian geodesic 
flow on the surface in this case it is necessary to determine de existence or non- 
existence of the solution of the equation (]45l) or (]4ll . We illustrate this case 
for the third-order surface 

(53) f(x)=xyz = c, c 7^ 0. 
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First we determine the dependence x = x(r, f,g), y = y(r, f,g),z — z(r, f, g). 
By considering that in this case 

g = { X yf + (xz) 2 + {yzf 

thus the functions /, g and r 2 are independent. Indeed if we introduce the 
cubic polynomial in Z : 

P(z) = Z 3 - r 2 Z 2 + gZ -f 2 = (Z- x 2 ){Z - y 2 ){Z - z 2 ), 

then by using Cardano's formula we obtain the dependence (|52|) . 

This case was examined already by Riemann in his study of motion of a ho- 
mogeneous liquid ellipsoid. More exactly, Riemann examined the integrability 
of the geodesic flow on (j551) . 

In [15] the author raised the problem. 

"Is it true that the geodesic flow on a generic third-order algebraic sur- 
face is not integrable?. In particular I do not know a rigorous proof of non- 
integr ability for the surface (|53|) " 

To prove the integrability of Lagrangian geodesic flow it is necessary to solve 
the non-lineal partial differential equation 

Gf(f, <£)(#||$x|| 2 - {yz** + xy^ z ) 2 ) = 2h(f)g 

Now we are not able to provide an answer to this question. 

6. DECARTES APPROACH FOR NON-HOLONOMIC SYSTEM WITH FOUR AND 
FIVE DEGREE OF FREEDOM AND TWO CONSTRAINTS 

In this section we apply Theorem 11.11 to study the non-holonomic system 
study in [8] and the well known non-holonomic system-the rattleback . 

Gantmacher's system 

Two material points Mi, M2 with equal mass are linked by a metal rod with 
fixed long and small mass. The system can move only in the vertical plane 
and so that the speed of the midpoint of the rod is directed along the rod. It 
is necessary to determine the trajectories of material points Mi, M2. 

Let (xi, y\) and (x 2 , 2/2) are the coordinates of the points Mi, M 2 . 
Introducing the following change of coordinates 

2ux = x 2 -xi, 2u 2 = yi-y2, 2u 3 = y 2 + yi, 2u 4 = x 1 +x 2 

we obtain the mechanical system with configuration space Q = M 4 , and La- 

1 4 

grangian function L = - u 2 — gu 3 . 

i=i 

The equations of the constraints can be rewritten as 



uiiii + u 2 u 2 = 0, U1U3 — u 2 Ui = 0. 



CARTESIAN APPROACH FOR CONSTRAINED MECHANICAL SYSTEMS . 23 

To construct Cartesian approach in this case we firstly determine the 1-forms 
Qj for j — 1, 2, 3, 4 as follow 

f2i = U\dui + u 2 du 2 , ^2 = u\du 3 — u 2 du 4 , 

Vt 3 = U\du 2 — U 2 dUi, fi 4 = U 2 du 3 + UidU4 

hence we obtain that T = u\ + u\. 

After some calculations we obtain that the vector field (T5]) takes the form 

v = v 3 (u x d 2 - u 2 d 1 ) + v 4 {u 2 d 3 + ui<9 4 ), T = (u\ + w 2 .) 2 , 

where i/j = A 3 (wf + u 2 .) for j = 3, 4. 

The 1-form associated to vector v is the following 

° — u 2 du% + Uidu 2 ) + v 4 [u 2 du 3 + u\du±). 

Thus the 1-form i v da admits the representation 

t v da = AifJi + A 2 tt 2 + A 3 tt 3 + A 4 fi 4 

= (ul + ul){- u 2 d U2 (^±^) - .A^) - 2^ 2 - vi)n x 

2 2 2 2 

+K + ul) (u 2 d U4 (^±^) - Ul d u /-±±^) + is 3 is,)n 2 

+(u{ + ul) Ui(u 2 d U3 v 3 + Mi<9 n4 z/ 3 + u 2 d Ul v 4 - uid U2 v^VL 3 

+(ul + ul)v 3 (u 2 d U3 P 3 + u x d Uik v 3 + u 2 d Ul p 4 - u^uAq^ 

If A 3 = A 4 = then u 2 d U3 v 3 + u x d Ui v 3 + u 2 d Ul u 4 - u x d U2 v 4 = 0. 

Cartesian approach generate the following differential equations respectively 

(54) -ui = -u 3 u 2 , ii 2 = v 3 u u ii 3 = u 4 u 2 , u A = v 4 u x 
and 

(55) u 2 d uz v 3 + u x d UA v 3 + u 2 d Ul u 4 - u x d U2 v 4 = 
It is easy to show that the functions v 3 , v 4 : 

, < , o 2\ /^( — 5^3 h) 9/9 9\ 

(56) ^ = <?3(«l+f 2 ). Ui= \ ( 2 i 2\ ^K + ^J, 

where 51, h are constants, are solutions of ( |55l) . f|54|) as a consequence 

2||v|| 2 = (u\ + u 2 2 ){v 2 3 + ^) = 2(-gu 3 + h). 
Under these restrictions Lagrangian approach generate the differential system 

(57) ill — Ai«i, u 2 = Aiu 2 , u 3 = —g + A 2 u 1 , u 4 = A 2 u 2 . 
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The solutions of ( 1551) are 

u\ —r cos a, u 2 = rsma, a = ao + g 3 (r)t, 
n 9 9 V^gC 

u 3 = ui H -— t ^— — sin 2a — — cos a, 

' ' 2g 3 (r) Agi(r) g 3 (r) 



u A = -h + + (-^- sin a + C)\ 

2g ^V2g 3 {r) 

where C, r, a > wjj, h, are arbitrary constants, #3 is an arbitrary on r function. 

To compare these solutions with the solutions obtained from the classical 
approach, we determine the equations of motion obtained from the d'Alembert- 
lagrange principle 

(58) Ui = UiUx, U 2 = UiU 2 , U 3 = -g + U4 = ~/i2«2 

where fxi, \i% are the Lagrangian multipliers. 

After the integration of the system (!58l) we obtain [8] 

/ / 

(59) ill = -tpii2, u 2 = <fUi, u 3 = -u 2 , u 4 = -Mi 

where (tp, r) are the polar coordinates: U\= r cos ip, u 2 = r sin <p and / is 
a solution of the equation 

(60) / = -^-u 2 

r 

The solution of (160|) is / = ^ C ° S — + 27 where 7 is an arbitrary constants. 

f 

Clearly if we choose v 3 = ip, v± = — then we the vector field associated to 

r 

system f)59p can be obtained from Cartesian vector field. 
The rattleback. 

The rattleback's amazing mechanical behaviour is a convex asymmetric rigid 
body rolling without sliding on a horizontal plane. It is known for its ability to 
spin in one direction and to resist spinning in the opposite direction for some 
parameters values, and for others values to exhibit multiple reversals. Basic 
references on the rattleback are [281 El E2l [2] . 

Introduce the Euler angles ip, (p, 9 using the principal axis body frame rela- 
tive to an inertial reference frame. These angles together with two horizontal 
coordinates x, y of the center of mass are coordinates in the configuration 
space Q = SO (3) x R 2 of the rattleback. 
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The Lagrangian of the rattleback is computed to be 
L — \(h cos 2 ip + I 2 sin 2 ip + m(T l cos 9 - ( sin 9f)9 
\(h sin 2 ip + I 2 cos 2 ip) sin 2 9) + I 3 cos 2 9)<p 2 
+\{h + mTl sin 2 9)ifj 2 + f (x 2 + y 2 ) 

+m(r 1 cos 6* — C sin #)T 2 sin + (Ii — I 2 ) sin 6* sin ip cos -000 
C cos 0<j>ip + m5-(ri sin 9 + C cos 6>) 

where h, I 2 , h are the principal moments of inertia of the body, m is the total 
mass of the body, 

Yi = £ sin ip + j] cos ip, T 2 = £ cos ip — r] sin -0 

and (£ = ^(9,ip), rj = r](9,ip), ( = ((9,ip)) are the coordinates of the point of 
contact relative to the body frame. The shape of the body is encoded by the 
functions £, rj and (. 

The constraints are 

x - ol x 9 - a 2 ip - a 3 <p = 0, y + ^9 + f3 2 ip + f3 3 (p = 0, 

where 

cti = (— ri sin — ( cos #) sin (p, a 2 = T 2 cos 6 sin + I\ cos 0, 

(?Q!l C?Ct 2 <9ct3 

a 3 = r 2 sin0+ (Ticostf - Csin0)cos0, pi x = — — , (3 2 = — — , /5 3 = -7—. 

<70 ocp ocp 

Clearly that the rattleback equations of motion in this particular case for- 
mally contain the equations of the heavy rigid body in the singular case 
m — >■ 0, mg — > I, I 7^ 

To determine Cartesian approach for the rattleback we first determine the 
1-forms flj, for j — 1, . . . , 5. In this case we determine as follows 

Qi = dx — a\d9 — a 2 dip — a 3 dcp, Q 2 = dy + (3\d9 + j3 2 dip + (3 3 dcp, 

Q 3 = d9, f2 4 = dip, Q 5 = d<p 

Hence T = 1 and the vector field v : 

(61) v = A3X3 + A4X4 + A 5 X 5 , 

where 

X 3 = a!d x - fady + de, X A = a 2 d x - /3 2 d y + d^, X 5 = a 3 d x - /3 3 d y + 8$. 

We now proceed to the consideration of the particular case for which £, rj and 
( admits the development 

£ = £o + e£i(0, V), V = Vo + em{9,i/}), C = Co + VO 
where £o, Vo, Co? are constants and e is a small parameter. Under this con- 
sideration we obtain that the Lagrangian function can be represented as follow 

L = L + eLi + e 2 L 2 . 
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Below we study the case when e = 0. 

Let (x 1 , x 2 , x 3 , x 4 , x 5 ) be a new set of variables derived from x, , y, 9, ip, <fi 
by the transformation 

ip = x 1 , <p = x 2 , 6 = x 3 , 
y + Co si n $ cos + cos 9 sin <p — T° sin <p = x 4 , 
x + Co sin 6* sin — cos 9 cos + T° cos (j) = x 5 , 

where = £ sin ip + r] cos ip, T° = £o cos ip — r] sin ip. 

The vector field v and the constraints on account of this change, take re- 
spectively the form respectively 

v = (a, b, c, 0, 0) x 4 = 0, x 5 = 

where a = a(x 1 , .., x 5 )), 6 = ^(x 1 , .., x 5 )), c = c(x 1 , .., x 5 )) are the C 1 functions. 

In the coordinates x = (x 1 , 5 ) the Lagrangian function L be- 

comes to the function 



L = i G jk x j x k + mgiT^smx 3 + Co cos x 3 ), 
j,k=i 

where G = (Gjk(x)) = (Gjk) is the Riemann metric. 

We shall now determine Cartesian approach under the given conditions. 

Proposition 6.1. The vector field v(x) = (a, b, c) is a Kummer vector field. 

Proof. Indeed the 1-form associated to the vector field v is 

a = p-^dx 1 + p 2 dx 2 + p 3 dx 3 , p k = G kl a + G k2 b + G k3 c, k = 1,2, ..,5 

then 

5 

i v da = Ajdx^ 

where 

dpi dp 2 ^, , ( dpi dp 3 dp 2 dp 3 dp 2 dp± 

A _ r a P3 d P2)h+( dp3_d Pl -_ d P\ _ d P* h _ d P* 



9pi <9p 2 , <9p 3 
A5 = -^ a -^ 6 -^ C 
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, x 2 


= b, x 3 


= c, 




dpi 




dpi 


dp 3 . 


dx 2 


dx 3 


dx 1>C 


dp 2 


dx 2 


,dp 2 


9pi )a 
dx 2 


dx 3 


W 


dp 3 


9p2 \ 7 . 

'dx 3)b + 


( dp 3 


dp x 
&> 


dx 2 




> C5), 


b = b(x\. 


X 2 X 3 


, C4, C 5 ), 



We have therefore that the differential equations generated by Cartesian ap- 
proach are respectively 

x 1 - 

Ai = (^-^f)6+(^-^)c = 0, 

(62) 

A 2 = ( _^__i_i )r+( _^_ LLl )a . (). 

where a = a(x 1 , x 2 , x 3 
Let rotv(x) be the vector field 

ro tv - 1 / dp 3 _ <9p2 <9pi_ _ 9^3 9^2 _ dpi T 
^JdetG dx 2 dx 3 ' dx 3 dx 1 ' dx 1 dx 2 
then the last three equations in fl62l) can be rewritten as 

[v(z) x rotv(x)] = 0. 
Thus the vector field v(x) = (a, b, c) is a Kummer vector field. □ 

For the general case, i.e., when the £, 77 and £ are functions on the variables 
9 and if) Cartesian approach produce the following equations respectively 

x = v(x) 



3=1 

E/up 3 upj ^ a ( UVi up i\ f3 L ( UVb UFj )) v "i=Q 
\ r>nr3 rinfl rlnrJ F)nr^ F)nf3 F)t^ ) 

where if) = x 1 , 
formula flBTl). 



(il! - ill _)_ a 2 (— - — ) - f3'A - — ) )v j = 

^ ^ \ rJ'rJ r)i™2 /-rVJ f)nr^ ft<r3 r)i"5 / 



dpi 
dx 3 


dpj 
dx 1 


+ a 2 ( 


dpi 
dx 3 


dp, j 
dx 4> 




<9p 5 
dx 3 


dpj 
dx 5 


dp 2 


d Pj 


+ oc 3 { 


dp a 


dp ^ 




dp 5 


dpj 


dx 3 


dx 2 


dx 3 


dx 4> 


dx 3 


dx 5 


'dp 3 _ 
.dx 3 


d Pj 
dx 1 


+ Cti 


dp 4 
[ dx 3 


d Pj 

dx A> 




,dp 5 
y dx 3 


d Pj 
dx 5 


<f> = 




9 = 




y = x 4 , 


a; 


= X 5 


and ' 



7. Inverse problem of dynamics 

Introduction 

This section is devoted to apply Corollary II .4l to study the problem of finding 
the field of force that generates a given (N — l)-parametric family of orbits 
for a mechanical system with N degrees of freedom. This problem is usually 
referred to as the inverse problem of dynamics. We study this problem in 
relation to the problems of Celestial Mechanics. 
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One of the fundamental classical problems in celestial mechanics is to de- 
termine the potential-energy function U such that every curve from a given 
family of curves will be a possible trajectory of a particle moving under the 

dU 

action of potential forces F, admitting U; i. e. F = — — . 

ox 

In the modern scientific literature the importance of this problem was al- 
ready acknowledged by Szebehely [3], [27] 

The first inverse problem in Celestial Mechanics was stated and solved by 
Newton (1687) and concerns the determination of the potential field of force 
that ensures the planetary motion in accordance to the observed properties, 
namely to Kepler's laws. 

Bertrand (1877) pQ proved that the expression for Newton's force of attrac- 
tion can be obtained directly from the Kepler first law to within a constant 
multiplier. 

Bertrand stated also a more general problem of determining a positional 
force, under which a particle describes a conic section under any initial condi- 
tions. Bertarnd's ideas were developed by [5] [26], [10], [6], and [7J. 

Dainelli in [5] essentially states a more general problem of how to determine 
the most general field of force (the force being supposed to depend only on the 
position of the particle on which it acts) under which a given family of planar 
curves is a family of orbits of a particle. 

The solution proposed by Dainelli is the following . 

Theorem 7.1. The most general field of force F— (F x , F y ) capable of gener- 
ating the family of planar orbits f(x,y) = const can be determine as follows 

0, m 



(63) F x = -\ 2 {f, d y f} - \{f, \}d y f, F y = \ 2 {f, d x f} + A{/, \}d x f, 

where {f, } = d x fd y — d y fd x and A is an arbitrary function which depends 
on the velocity with which the given orbits are described. 

By considering that the components F x and F y are to be functions of the 
position of the particle, we can take A to be an arbitrary function on x and y. 

The above expressions for the field of force under which the curves of the 
given family are orbits were first given by Dainelli [5]. 

After some calculations we can prove that (1631) can be rewritten as follows 

<9-||v|| 2 

(64) F = - \ (d X (Xd X f) + dy(\dyf)) ^ 

where v = (-Xd y f, \d x f) . 

Suslov in [26] stated and solved a problem which was a further develop- 
ment of Bertrand's problem. He shows that, given a (N — l)-parametric 
family of orbits fj = fj{x) = Cj for j = 1,2, ... N — 1 in the configuration 
space of a holonomic system with N degrees of freedom and a kinetic energy 
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1 N -1 

T = — Gjk{x)x 3 ' x k = -| |x| | 2 , it is necessary to determine the potential field 

j,k=i 

of force under which any trajectory of the family can be traced by the repre- 
sentative point of the system. Suslov deduced the following system of linear 
partial differential equations with respect to the require potential function: 

89 dU 39 dU U + hf 39 39 39 39 



3A k 3x N 3A N 3x k 9 \3A N 3x k 3A k 3x N 

+ V A m ( 99 dH - " dH ) 

^ { 3A k 3A N 3x m 3A N 3A k 3x m> 

for k — 1, 2, .., N — 1. where 9, A 1 , A 2 , A^ are functions: 



k=l k,j=l 
N 

A k = J2 G Ax)& j , k = 1, 2, ., N, a = 1, 2, ., N - 1, 

3=1 

and proved that theses equations represented the necessary and sufficient 
conditions under which the equations of motion of the study mechanical system 
admits the given iV — 1 partial integrals. 

Assuming that given trajectories admit a family of the orthogonal surfaces, 
Joukovski in [10] constructed the potential-energy functions in explicit forms 
for systems with two and three degrees of freedom. 

The following theorem was enunciated by Joukovsky in 1890 

Theorem 7.2. If q = const is the equation of the family of curves on a 
surface, and p = const denotes the family of curves orthoqonal to these, then 
the curves q = const can be freely described by a particle under the influence 
of forces derived from the potential- energy function 



V = A 1 (p)(g(p) + J Hq)§- q (^)dq 



where h and g are arbitrary functions, and Aj denotes the first differential 
parameter. 

A new approach to the problem of constructing the potential field of force 
was proposed by Ermakov in [6], who integrated the equations for the potential- 
energy function for several particular cases. 

In the most general form the inverse problem in dynamics was studied in [2"H 
[22] . By applying the results presented in that work we propose the following 
new results: 

(i) Statement and solution of inverse Dainelli's problem for a mechanical 
system with iV degree of freedom. 
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(ii) New approach to solve the Suslov problem. 

(iii) Statement and solution of inverse Joukovski's problem for mechanical 
system with N > 3 degree of freedom. 

(iv) Generalization of Theorem 17.21 for mechanical system with N > 3 de- 
gree of freedom 

(v) Statement and solution of inverse Stackel's problem. 

(vi) General solution of Bertrand's inverse problem. 

The results listed above are obtained by applying Corollary 11.41 

Statement of the generalized inverse Dainelli problem. 

Given a iV — 1 -parametric family of orbits fj = fj{x) = Cj for j = 
1,2, . . . , N — 1 in the configuration space Q of a holonomic system with N 

degrees of freedom and kinetic energy T = - 1 |x| | 2 . Generalized Dainelli 's prob- 

lem is the problem of determining the most general field of force that depends 
only on the position of the system under which any trajectory of the family 
can be traced by a representative point of the system. 

Solution of the generalized inverse Dainelli problem 

The following proposition provides a solution to the problem above 

Proposition 7.3. Given a mechanical system with a configuration space Q and 

a kinetic energy T = -||cc|| 2 ; then the most general field of force F that depends 

only on the position of the system and is capable of generating the given orbits 
fj{x) = Cj, for j = 1, . . . , N — 1 where fi, ■ ■ ■ , /jv-i are independent functions 
can be determine from the formula 

d[\\\v\\ 2 ) N-i „ 

where 

(66) v = -\ N l(hlll±p=h*L =\{f l ,..., f N _ u *}, 

1/1) • • • ) JN-I, jn\ 
a^j = a N j(x) for j = 1,2, . . . N are functions: 

a Nj = (-l) N+j - l da A dfx A df 2 A ... A df^ x A df j+1 A ... A df N ^{d u ...,d N ), 

and \n and /at are arbitrary functions such that {f\, . . . , fN-i, /at} 0. 

Proof. Let us suppose that are given the N—l parametric family of trajectories, 
hence (d x fj, x) = for j = 1, ... iV — 1. In view of independence of functions 
fi, . . . , /jv-i we can solve these equations respect to velocity, thus we obtain 
the system x = v(x), where v is determine by the formula fl6"6"|) . 

After covariant derivation we obtain the equations of motion of the mechan- 
ical system (see proof of Theorem 11.11) 
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I -llvl 



ddT dT _ " V2"'" J x sp a d fi _ F 
rft 9x <9x <9x N ^ <9x ' 

therefore the proposition has been proved. □ 

The following proposition shows that Theorem 17.11 is a particular case of 
Proposition 17.31 

Corollary 7.4. For N = 2 and Q = M 2 the force field (|65|) coincides with the 
solution proposed by Dainelli . 

Proof. Indeed for N = 2 the field of force F takes the form 

ax ox 
On the other hand by considering that v = {—Xd y f, Xd x f) thus 

|V|| 2 = (XdJ) 2 + (Xdyf) 2 , 02! = OV^, dy) = d X (XdJ) + dy(Xdyf) 

hence we obtain the formula (l64l) . □ 
Corollary 7.5. For N = 3 the force field (1651) tofces £/ie /orm 

^(-ll^ll 2 

(67) F = K2 dx 7 + A (dh(rotv)df 2 - df 2 {rotv)df x ) 
Example 

Given a particle with Q = R 3 and kinetic energy T = -(£ 2 + if + £ 2 ). 
Construct the field of force capable of generating the two-parametric family of 

trajectories defined as intersections of the two families of surfaces 

(68) /i = C = ci, f 2 = H(t, V ,() = c 2 . 



The solution of this problem can easily be derived from Corollary 17.51 The 
vector field v, rotv are the following 



A 



dH d dH d 



drj <9£ <9£ dn 
d f,dH\ d 6 ( \dH\ d d 

rotv = d([ X ^)^ + dc{ X ^)^'^dC 
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Hence the require field of force is such that 



(69) fl » ^ 



c??7 \ 2 \ d£ drj J J drj 

C = o, 



^ d^\d^) drj \ drj ) 
In the next section we make use of the solution of the generalized Dainelli 

inverse problem for solving the Suslov and generalized Joukovski problems. 
Statement of Suslov 's problem 



Given a N — 1 -parametric family of orbits fj = fj(x) = Cj for j = 
1,2, . . . , N — 1 in the configuration space Q of a holonomic system with N 

degrees of freedom and kinetic energy T = -||x|| 2 . Suslov's problem is the 

problem of determining the potential field of force that under which any tra- 
jectory of the family can be traced by a representative point of the system. 

Solution of Suslov's problem 

We now propose a new solution to the Suslov problem. This solution we have 
obtained as a special case of the previous solution to the generalized inverse 
Dainelli problem. 

Proposition 7.6. Given a mechanical system with a configuration space Q 
and a kinetic energy T = -||a;|| 2 , then the potential field of force F = , 

capable of generating the given orbits fj(x) = Cj, for j = 1, . . . , N — 1 , can be 
determine from the formula 



d I -\\v 



dU " V 2 """ J x d f: 



dx dx ^ ^ ^ aNj dx ' 

j=i 



13 



N-l 



if and only if 

(71) A ^2 a N j(x)dfj = dh(f 1 ,f 2 , f N -i) 

3=1 

Clearly if ( ITT]) holds then the potential function U is such that 

U(x)= l -\\v\\ 2 + h{f 1 ,f 2 ,...,f N _ 1 ), 
where v is determined by the formula (|66|) . and X is an arbitrary function. 
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au N ^ 

Proof. From fl65|) follows that F = — — , if and only if A aNj(x)dfj is exact 

1-form dh. In view of the relations dfj(v) = we deduce that dh(v) = 0, 
hence in view of independence of functions fj for j = 1, 2, . . . , iV — 1 we get 
h= (/i,/2,...,/jv-i)- 

□ 

i2rampZe(Ermakov's problem). 

Given a mechanical system with configuration space Q = R 4 and kinetic 
energy 

T = Q(miOi + Vi) + rn 2 (x 2 2 + y^j . 

Construct the potential field of force capable of generating the three-parametric 
family of trajectories defined as intersections of the families of hyper- surfaces 

h = %1+Vi = ci, h = x l + vl = c 2, h = (x 1 -x 2 ) 2 + (yi-y2) 2 = c 3 

Corollary 7.7. Under the assumptions ofProposition^TTBforN = 2 we obtain 
that 

< 72) i^ + MWl + s.w)! 

if and only if 



(73) A (d X (Xdj) + dy(Xdyf)) df = dh(f) 

where v = (-Xd y f, \d x f) 

Another interesting application of the solution to the generalized Dainelli 
problem is the determination of the solution of the generalized Joukovski prob- 
lem. 

Statement of generalized Joukovski problem 

Given a N — 1 -parametric family of orbits fj = fj{x) = Cj for j = 
1,2, . . . , N — 1 in the configuration space Q of a holonomic system with iV 

degrees of freedom and kinetic energy T = -||x|| 2 . Assuming complemen- 

tary that the given trajectories admit a family of orthogonal hyper-surface 
5* = S(x) = cn then Generalized Joukovski problem is the problem of deter- 
mining the potential field of force that under which any trajectory of the family 
can be traced by a representative point of the system. 



Solution of Generalized Joukovski problem 
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Proposition 7.8. Given a mechanical system with a configuration space Q 
and kinetic energy T = -||i|| 2 , then the potential field of force F = , 

capable of generating the given orbits fj(x) = Cj, for j = 1, . . . , N — 1 , which 
admit a family of orthogonal hyper-surface S = S(x) = cjy can be determine 
from the formula 

dx dx dx ' dx dx dx dx 

if and only if 

( 74 ) jr^ dS ~ HttII v = -2dh(h, f 2 ,..., /*_x) 

ox ox ox 

where v is an arbitrary function. Clearly if ( 1741) holds then the potential 
function U can be determined as follow 



C7 =l,Vf ll a5 111 +ft(/ " h -' h -' 



( dS df 

the given trajectories then the following relations hold I — — , — — ] = for 

V ax ax 



Proof. In view of the condition that there is an orthogonal hyper-surface to 
the given trajectories tht 
j = 1,2, ...,JV- 1, thus 

(99 

(75) P G- 1 — = {f 1 ,f 2 ,...J N - 1 ,*} 

where p is an arbitrary nonzero function and G^ 1 = (G jfc ) is the inverse matrix 
of the Riemann metric G. Hence we obtain that the vector field fl66l) takes the 
form 

v = vG — , 
ax 

where v = Xp, therefore the 1-form associated to this vector field is such that 

fdS , \ 
a = v \ — — , ax consequently 
\dx J 

dv dS OS 
t v da = dv{y)dS — dS{y)dv = ^(7^—, —)dS — i>| — — \\ 2 du 

ax ax ax 

hence we easily obtain the proof of the proposition (see for more details proof 
of Theorem (II .ip ) . □ 

As a first application of the above proposition we have the following results 

d$(S) 

Corollary 7.9. If in Proposition 17.51 we suppose that v = then the 

dS 

potential function U can be determine as follows 

where ho is an arbitrary constant. 
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Proof. Indeed if v = — — — then then the 1-form a associated to vector field 

dS 

v is exact, thus i v da == 0, consequently (see formula flUj) ) dh = 0. Hence 

□ 

Example{ Bertrand's Problema [1]) 

Given a particle with configuration space Q = M 3 and kinetic energy 

Construct the potential field of force capable of generating the two-parametric 
family of trajectories defined as intersections of the two families of surfaces 

fi = C = c 1 , f 2 = Ve + tf + b£ = C2. 
From (175|) we obtain that the function S : 

dS_ _ 1 r] _ 1 / g _ 

d£ ~ Pvf+?' Sly ~ p I'v^V^ 



Hence, by choosing p = r] we obtain that S 1 = ln(£ + a/^ 2 + rf) — (b + 1) In r] 
and ^ = rj\. Clearly the field of force is potential in particular if v = $(S). 
The general solution of this problem we give below. 

Example 

Given a particle with configuration space Q = R 3 and kinetic energy 

Construct the potential field of force capable of generating the two-parametric 
family of trajectories defined as intersections of the two families of surfaces 

fl= XZ = Ci, f 2 =yz = C 2 

From (1751) we obtain that the function S = -(x 2 + y 2 — z 2 ), thus the condition 
( I74p takes the form 

(xd x v 2 + ydyU 2 — zd z u 2 )(xdx + ydy — zdz) — (x 2 + y 2 + z 2 )dv 2 

= -2dh(f 1 ,f 2 ). 

There are two obvious solutions v = v[x 2 + y 2 — z 2 ) and v = z. The first 
solution produces the potential function U = U(x 2 + y 2 + z 2 ) — ho, h — h 
which coincide with the solution obtained by Joukovski and the second gives 

the potential U = -z A — ho, h— -(f 2 + /f) — h . 
The following result is a generalization of Theorem 17.21 



36 CARTESIAN APPROACH FOR CONSTRAINED MECHANICAL SYSTEMS . 

Corollary 7.10. If x 3 = Cj = const, for j = 1,2, ..,N — 1 are the equations 
of the N — 1 parametric family of curves on Q, and x N = const denotes the 
family of curves orthogonal to these, then the curves x J = Cj = const can 
be freely described by a particle under the influence of forces derived from the 
potential- energy function 



U 



1 ( g{x N )+ srj h{x\x\..,x N - 1 ) dGNN{xl dx f 



Gj\[j\f (x^~ , x 2 , ..,x j i 



where h = h(x l , x 2 , . . . ,x N *) and g = g(x N ) are arbitrary functions. 

Clearly, forN = 2 we obtain the Joukovski theorem given in the introduction. 

Proof. Whereas in the study case fj = x\ j = 1, 2, . . . iV — 1 and S = x N then 
from Proposition 17.81 we obtain that the relations hold 

" = , ..O....Y ,. 

OX 3 OX 1 

d2U _ d (y 2 G NN ) 




dx N dx N 

where Gkj = Gkj{x x , x 2 , . . . x ) for k,j = 1, 2, . . . N, if and only if 



X - .2 



(76) (J2 G jN ^j)dx N - G NN dv 2 = 2dh(f u f 2 , . . . , J^), 



dx j 



assuming that the Riemann metric is orthogonal then 

/ ; . .i / < \ \ ^ 



dv 2 = -2G NN dh, where 



G 



NN 



v 2 = 2g{x N ) -2 J G NN dh 

N-l 

^+E 2 / 

3=1 ^ 



2^) - 2G^ + y> / h ^r dxj ' 



where g = g(x ) is an arbitrary function. 

1 i/ 2 
Consequently in view of the formula J7 = -v 2 G NN + h = — — h ft we 

2 2Gnn 
obtain the proof of the proposition. □ 

Now we apply Theorem 17.81 to solve the inverse problem which we will call 
the inverse Stackel problem. 

Statement of inverse Stackel problem. 

Given aiV-1 -parametric family of orbits 

(77) L = Ux) = V / lf / k ^ X % dx k = c», n= 1,2,..., N- 1, 
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where K k (x k ) = 2^ fc (x fc ) + 2^ otj(fkj(x k ), a>k, for k = 1, 2, ..iV are constants, 

3=1 

in the configuration space Q of a holonomic system with iV degrees of freedom 
and kinetic energy 



(78) 



1 = - X 

2 M I 



1 - 



2 ^ A? 



where A- 7 = for j = 1, 2, . . . , iV are functions such that 



iV 



A' 



d<£ a = y^V?fcq(^)^ fc , V'fca = (fka{x k ), for = 1, TV, a = 1, . . . , N are 

k=l 

arbitrary functions. 

The inverse Stackel problem is the problem of determining the potential 
field of force that under which any trajectory of the family can be traced by a 
representative point of the system. 

Solution of inverse Stackel problem 



Proposition 7.11. Given a mechanical system with a configuration space Q 

dU 

and kinetic energy (178jl . then the potential field of force F = ~q~- capable of 
generating the given orbits ( 177]) is the field of force with potential function 



(79) 



U = v (S) [ r- + a x ) - he 



r N N 

where S = 2 / \^ . ^ k (x k ) + atjipkj{x k ) dx 

J 3=1 \ 3=1 



Proof. In view of the equality 



{fi, h, for-h *} 

{fl, f2, fN-1, /w} 



q 1 {x 1 )d<4) 1 {d 1 ) 



q N {x N )d(p l {d N ) 



q N (x N )d(p N ^(d N ) 







N 



N 

E 



A> 



N 



3=1 



dxi 



x k ) 
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where qj{x^) = — , consequently from (175|) and fITTj) follows that 

\/K k (x k ) 

1 dS 1 dS 
v = -\ N G — = vG — , 
ox ox 

hence 

N N / N 

||v|| 2 = ,^.V'(/x ; iV)) 2 = v 2 ^A k h® k (x k ) + 2^> jW 

i=i ^=1 V j=i 

TV AT TV 

= 2v 2 A k Mx k ) + 2^ 2 a i E ^W**) 

k=l j=l k=l 

= 2z/2 I {V?l, V?2, • • • , V?7V-1, ^} ^ {<?!, V?2, <PN-1, Pj} 

\{<Pl, <P2, ■ ■ ■ , <PN-1, <Pn} ~[ j Wl, <P2, ■ ■ ■ , fN-1, Pn} 

= 2u 2 ( ill' <P2, ■ ■ ■ ; ^} + a \ 

\{<Pl, ¥>2, • • • , <Pn-i, ¥n} V ' 
AT 

where eM> = ^ ^ k (x k )dx k . 

3=1 

On the other hand if we choose v = u(S) then from Corollary 17.91 follows 
the proof of the proposition. □ 

Corollary 7.12. If in the previous proposition we suppose that v(S) = 1 and 
a>i = ho then the potential function (|79p coincide with the Stackel potential [I] . 

Proof. Indeed, if v(S) = 1 and ot\ = h then [I] 

u = M.y>2,...,^-i,^} = y A k Mx% 

□ 



Statement of Bertrand's problem 

Let a particle with configuration space Q = R 2 and kinetic energy T = 

1 

— (x 2 + y 2 ) . Bertrand's inverse problem is the problem of determining the 

most general potential field of force capable of generating the one-parametric 
family of conies / = \J x 2 + y 2 + bx = c, where b is the eccentricity. 

Solution of Bertrand's inverse problem 

Now we study the problem of constructing the potential field of force, which 
is capable to generate given conies. We prove the following proposition. 
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Proposition 7.13. The potential- energy function U capable of generating a 
one-parameter family of conies with eccentricity b is the function 

U = a_i#_i(cos 0) + logr(l + bcos8) 
iez\{-i} v 

if b 7^ where a,j j G Z, ifj are real constants and Hj, j G Z are solutions of 
the Heun equations with singularities at the points 



1 + 6 
0,1,^,00 



and with the exponents 



(0, 26 ) ; t '- 7 2fe ^' ^ ) ^ + l >^~ 1 ~^ 1 ~ 3 > 

and 

if b = 0, where \& = ^(cos^) ana 1 /i = /i(r) are arbitrary functions. 

Proof. From Proposition 17.71 follows that the require potential field of force 
exist if and only if 

x \ dX 2 y dX 2 2X 2 dh 
+ b\—+ , = — + = 2— 



a/x 2 + y 2 J dx ^/x 2 + y 2 dy r df 

By introducing the polar coordinates x = rcos9, y = rsin(9, we find that 
previous condition condition takes the form 

(1 + bcos6) — - bsinddX2 | 2X2 - 2 dk 
dr r 89 r df 

or, equivalently, 

(so, (1+6T) f! + M^i!)|! + ^. 2 | 

or r or r oj 

where / = r(l + br), r = cos#. We embark now upon the study of the case 
when b ^ and h is such that 

(81) h(f) = v_ 1 \n\f\+Y,^TTT> 

where Uj, j G Z, are real constants, and A is determined in such a way that 

(82) A 2 = 5>i(r)lTi(T) 

It is clear that the series (15TD and (1521) are formal series. 
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By inserting ( IHTj) and ( 1821) into ( IHUl) we obtain 

£ + *)*^W> + *M ( 6(1 _ T ^ + 2 ) _ 2 ^ (1 + bT)1 ) 

jez ^ V / / 

This equation holds if 

if)j(r) = ajr J+1 , v, = -ajKj 
for j e Z and we determine iij as a solution to the equation 
(83) 6(1 - r 2 )#;(r) + ((j + l)br + j + 3)^(r) + 2Kj{l + 6r) J = 
for j G Z. 

The general solution of this equation is 

21^ r (1 + br)'J 
b J (l-r 2 )^(r)' 
3 + 1 j + 3 J + 1 J + 3 

Ci(r) = (l-r) 2 + 26 (r+1) 2 26 

where Cj, j e Z are arbitrary constants. 

Under these conditions, the required potential-energy function U results in 
the form 

U(r, t) = ^A 2 (l + b 2 + 2br) - h(f) = £ ajU^r, r) 

where 

U 3 (r,r) =l r i+ l H J (T)(l + b 2 + 2bT) + -^-P +1 ifj^-1 

tf-i(r, r) = ^-i(r)(l + 6 2 + 26r) + K_ x In |/|. 

We will study the subcase when 6=1 separately from the subcase when 6^1. 
If b= 1, 

tf(r, r) = A 2 (l + r) - h(f) = a Mr, r). 



where 



t) = r^(i - r)^+ 2 (q, - 2 Kj J { iXy+* dT ) + jTT /J+1 ' lf 3 * - 1 

tf-x(r,r) = (1 - r) (C7_! - 2^ / (1 _ r) t (1 + r) ) + ^ H/l- 

Easily verifies that 

C - 2 o g -V/' dT 1 ^ C7_, JT_ 2 
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1 - r 

where g(r) = In \ j - — — . Therefore, if 6 = 1, 



If 6 1, 6 0, it is easy to prove that 

1 - b 



H _ 2iT) = <izlL»^ 2K - 



1 + 6 ~ z (6t + 1)(1-6 2 ) 
(1 + r) 26 

C/-2 r, r = — — 1 + 6 2 + 26r - = + G ' 

2r r(or+l) r(tr — 1) r 



where 



1 



gw = ^(ttt)^ (1+62 + 26t) 

Under these conditions, the potential function U takes the form 
C/(r, r) = — — - G(t) + _ + ^ a 3 U 3 {r, r) 

Summarizing the above computations, we deduce that if 6 7^ the function 
[/ is represented as follows: 

tf(r,T) = £ + ^+5>tfi(r,r) 

jy-2 

where a is a constant and /3 = /?(t) is a certain function. 
If 6 = 0, then f = r and condition fISUl) takes the form 

\2 

<9 r A 2 + 2— = 2d f h(f) . 
r 

Therefore, 

r 2 A 2 = 2 J r 2 d r h(r)dr + 2^(T), 
which rearranged results in the expression: 

A 2 = 4- / r 2 d r h(r)dr + ^All = 2/i(r) - / h(r)rdr + ^ ' ' ' 



where \I/ is an arbitrary function. 
Hence, 



f/(r,r) = ^-- / 6(r)dr. 
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To end the proof of the claim will establish the relationship between the 
functions Hj for j e Z and the solutions of Heuns equations. 

The canonical form of Heun's general equation will be taken as [23J 



d 2 x 



cb 2 \2 2—1 z — ajdz z(z — l)(z — a) 

here x and 2 are regarded as complex variables and a, f3, 7, 8, e, a, b are param- 
eters, generally complex and arbitrary, with the only condition that a 7^ 0, 1. 
The first five parameters are linked by the relation a + ^ + 1 — j+8 + e. 

The equation is, therefore, of the Fuchsian type , with regular singularities at 
the points z — 0, 1, a, 00. The exponents at these singularities are, respectively, 
(0,l-7),(0,l-e),(0,l-<y),(a,/3) 



dx 



a(3z — B 



0, 



Now we establish the relation between equation (J83J) and Heun's equation. 
By fulfilling the replacement 

we can easily obtain the following representation for (1851) : 

z ( z - i)¥k + + b - 2bz)(j + 1) + 2)H j {z) - ^(1 + b - 2bz) j = 

dz 2b ' b 

By differentiating and fulfilling some straightforward calculations, we deduce 

that the functions 

z ( z -l)^+l((l+b-2bz)(j+l)+2)H j (z)yi+b-2bz)^ (j e Z) 

are first integral of the Heun equation: 



F 3 {z 



(85) 



d 2 Hj 
~Jz T 



(l 
\ 

Dz 



-a(l + j)-~ a(l+j) + ^-j 



z - 1 



z(z — l)(z — a) 



where a 



l + b 
2b 



z — a 



dIL, 

dz 



By comparison with classical Heuns equation, we obtain: 

1 + 6/ \ 1 P l + b, s 1 

T, = i-—a +} )— s, = j. 



B; 



J 



1 

1 + 6,. 



-(l + e ,)(2- 7j -^), 



e i = -3 



b + 2b 



(f - 1) = -o(2 - 7i - Sj) - (1 - 7ih 



Evidently, when the given conies are parabolas then in (185]) we have the 
confluence of singularities. In fact, in this case b = 1 so a = 1, and as a 
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consequence the Heun equation is transformed into hypergeometric differential 
equation 

d 2 Hj (-1-3 2 - A dHj (j 2 -1){z-1)-j TT . . n 

for j G Z. Which completed the proof of Proposition. □ 
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